NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



4k* A 


79-FM-31 


JSC- 16079 


Lyapunov Stability and its 
Application to Systeirs of 
Ordinary Differential Equations 



(5 IS A-’** -= *'?>.) LYAPUW3V STAPILITY 
APPLICATION 70 SYSTZ5S OF CPDTS1PT 
DIFFFrEy’IAL E^UATIOSS (7ASS) 32 p 
nC A 03,'* 7 AOT 


AUC IIS 
CSCL 1 2 A 


G 3/6 4 


730-13335 

"Tildas 

19006 


» 

» 


Mission Planning and Analysis Division 
September 1979 



/y 


A. p 5 ^ ^ 

NASA 

§ £ &£? A 

[to ' C ^ -o 

v*i -f 

National Aeronautics and 
Space Administration 

V 3 ^ 'T'V 

V vj r ; f 

Ntf /v- - • ? 


\ 


Lyndon B. Johnson Space Center 

Houston Texas 





79-FM-3' JSC-16079 

SHUTTLE nOGUH 

LYAPOTDf STABILITY AH) ITS APPL1CATIOB TC SYSTEMS OF 
ORDHABY D IF F MMTI AL BQUATICMS 

By Suel H. Kennedy* 

Technology Development Office 


Approved: 




/ 


Victor R. Bond 

Technology Development Office 


Approved 


, Chief 

Mission Planning and Analysis 




Mission Planning and Analysis Division 
National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, Texas 
September 1979 

•Summer Faculty Fellow, California Polytechnic State University at San Luis 
Cbispo 


79FK21 


A'~n»WLEDGMorr 

The author extends his sincere appreciation of the invaluable assistance of 
Yictor Bend of the Mission Planning and Analysis Division at Johnson Space 
Center. Of exceptional value were his numerous and informed discussions 
on the various formulations of the equations of orbital action. 


PRECEDING PAGE ELANK NOT FIL: , 5 


m OT F?L » , ' 




79FM31 


C0STEN7S 

Section Page 

1.0 INTRODUCTION 1 

2.0 LYAPUNOV STABILITY AND ITS GEOMETRY 2 

3.0 LINEAR SYSTEMS 8 

3.1 PERTURBED LINEAR SYSTEMS 15 

4.0 PERTURBED NONLINEAR SYSTEMS 17 

5.0 CONCLUSIONS 27 

6.0 REFERENCES 28 


v 


79FM31 


FIGURES 

Figure Page 

1 A Lyapunov tube 3 

2 A second trajectory within a Lyapunov tube 4 

3 Solutions of y = 1 - y 2 21 

4 Solution of u = 2u - u^ 21 

5 Solutior.3 of 6y = 26y 21 


vi 


79FM31 


1.0 INTRODUCTION 

Two of the more influential factors in the study of numerical integration of a 
system of differential equations are mathematical (Lyapunov) stability and numer- 
ical stability. Some of the confusion surrounding these concepts seems to re- 
sult from certain similarities; specifically, factors such as (1) both types 
of stability are related to the behavior of the eigenvalues of the system under 
analysis, (2) similar terminology, (3) interrelationships between the concepts, 
and (4) instability of either type frequently results in an unacceptable loss of 
accuracy of numerical approximations. 

Lyapunov stability is concerned with the qualitative behavior of the solutions 
of a vector differential system 


2 = I( t ,i) 


in some neighborhood of a specific solution 2o^ of the same system. Appro- 
priate behavior on some infinite interval t 0 < t < °° will result in 2.0^) 
being classified as stable, unstable, asymptotically stable, or one of several 
additional categories. Lyapunov stability is a global property of a particular 
solution (determined by the nature of its relationship to solutions in its imme- 
diate vicinity) . 

A useful technique in the study of the stability of a particular solution y_o(t) 
is to linearize the function f(t,2) about the solution ^(t); i.e., to trans- 
late the origin to y_ 0 (t). If the nonlinear perturbation is admissible in some 
sense, then the stability nature of the linear system can be transferred to the 
solution 2o^ • The eigenvalues of the linear system are characteristic for 
the classification of its solutions (with regard to Lyapunov stability) . 

This report provides an outline and a brief introduction to some of the con- 
cepts and implications of Lyapunov stability theory. Various aspects of the the- 
ory are illustrated by the inclusion of eight examples, including the Cartesian 
coordinate equations of the two-body problem, linear and nonlinear (Van der 
Pol's equation) oscillatory systems, and the linearized Kustaanheimo-Stiefel ele- 
ment equations for the unperturbed two-body problem. Lyapunov's direct method 
was not considered, nor was an attempt made to present the most relaxed version 
of a particular theorem. For additional information at the introductory level, 
see references 1 and 2. At the advanced level, references 3 and 4 provide an 
extensive survey of the field. 

It is anticipated that this report will provide contrast to a second report on 
numerical stability. 
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2.0 LYAPUNOV STABILITY AND ITS GEOMETRY 

Consider the concept of Lyapunov stability for arbitrary systems 


i = f (t,£) 


( 1 ) 


where 


Z = <y i C t.) y n (t)) T 

• _ d £ 

1 = — 
dt 


and 


f(t,z) = (fi (t,z) » • • • t 

It is assumed that f is defined and continue 's in some tube in t x R n space. 


T = ( ( t ,y) : £ t < -, Mzll < a 0 > 

where a 0 is a given positive constant, t 0 is an initial time, and | j j | de- 
notes the euclidean norm 


llzll =yjl vi* (2) 

or some equivalent norm such as 

I I zl I = raax t|yiU (3) 

i = 1 , . . . ,n 


The continuity condition on f ensures the existence of solutions to initial 
value problems associated with equation (1). In many cases, a Q = 00 for 
the tube T. 

Let Zo^) 2 Z^tJto’Xo^ denote some solution to the initial value problem 
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i = l 

y.(t Q ) = Zo ) (**) 

then 

Definition 1 . A solution ^o(t) is said to be stable in the sense of Lyapunov 
if for every e > 0 there exists a 6 > 0 such that if £•) is any solution of 
equation (1) with ^i(t 0 ) = Xl > anc * the number 

IIZl - Zoll <5 (SI) 


then 


Izi(t;t 0 ,j£i) - Z 0 ( t ; t o*Z 0 )|l < £ 


(S2) 


for all t, t ^ t 0 . If ^ Q (t) is not stable, then it is unstable. If ^ 0 (t) 
is stable, and if the additional condition 


lim | |Zl ( t» to»Zl ) - Zo(t;t 0 ,iT 0 ) | | =0 (S3) 

t-K= 


is valid, then &>(t) is asymptotically stable. 

Discussion . In general, stability or asymptotic stability is a property of a 
single solution io(t;t 0 ,£ 0 ) and requires a certain behavior over an infinite 
time interval (t 0 ,®). In geometric terms, consider a tube about the solu- 
tion Xo(t;totZo) °** ra di- us e > 0 (o arbitrary). 

Figure 1 represents a Lyapunov tube about a vector £ 0 ( t » t o>Zo) in t x R^. 
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If is staDle * then there exists a smaller neighborhood of radius 

6 of such that any solutions £i(t; t 0l £i) starting in that smaller 

neighborhood will never stray farther than e distance from Xo(t;t 0 ,£ 0 ) (i.e., 
solutions starting in the o disk never leave the e tube). (See figure 2.) 
Asymptotic stability merely ensures that all solutions starting in 
some 5 neighborhood of ^o tend to the solution y.o^t;t 0 ,^ 0 ) as t -*■ qo. 



Figure 2.- A second trajectory within a Lyapunov tube. 

The first example will establish that circular trajectories of the Cartesian co- 
ordinate system formulation of the two-body problem are unstable. 

Example 1 ■ For the two-body problem with masses m and M subject to the 

k 2 mM T 

universal law of gravitation — where r = C y i ,y2 *y 3^ 1 » r = |£.| (all 
functions of time t) , the resulting differential equation is 

r + — r = 0 ( El . 1 ) 

- r 3 " 

where 

p = k 2 (M + m) = gravitational parameter 
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As a first-order vector system, equation (E1.1) can be expressed as 



where 

v = r = (y-|(t) ,y2(t) ,y 3 ( t ) ) T 


03x3 is the 3 by 3 zero matrix, and 13x3 is the 3 by 3 identity matrix. Sys- 
tem equation (El. 2) is clearly nonlinear, and a well-posed initial value problem 
would require initial positions and initial velocities; i.e., 



where (£o»Xo^ are 6i ven * Because motion occurs in a plane, it can be assumed 
(without loss of generality) that y3 = 0, and the motion is in the (y-|,y 2) 
plane. Equation (E1.1) implies 


.. u 

y 1 + ~ y 1 = 0 

r-J 

U 

y 2 + y 2 = 0 
r - 5 


(El. 4) 


where 


r 


W 1 2 ♦ y2 2 


Next, change to polar coordinates by setting 
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= r cos <fc 
y 2 = r sin <c 

where r,4> are functions of t. This leads to the well-known and equivalent 
system 


r - r^ + - = 0 

r (El. 5) 

rtf + 2r<j> = 0 

For circular orbits, r is constant; hence, r = r = 0. Consider two orbits of 
fixed radii r-| and r2, then from equation (El. 5) 


y 




/y t 

1 ' 


i = 1,2 


Thus, two solutions of the vector formulation of equation (El. 4) are given by 


ni 

and 


' /y t 
cos 

r-j 3/2 


sin 


/£ t\ 

r~) 
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J 


Li = r 2 


cos 


/y t 


/y t 


r 2 


3/2 


sin 


r 2 


3/2 


Place r ] in Hi, but consider ]]2 to vary as a function of r 2 - In order to 
show instability of Hi , it is required to show that there is some e > 0 such 
that for every 6 > 0 there exist solutions for which 


1 1 D.1 ( t ; t 0 , r 1 ) - H 2 C t; t 0 ,r 2 ) 1 1 > e (El. 6) 

for arbitrarily large values of t, and | | r -j - r 2 1 S < That is, every 6 

disk about r-j contains solutions (starting in that disk) that escape the 
Lyapunov tube for arbitrary large t. An easy computation shows that 


| J! -j - 1 1 = *^2* | 1 - r]T 2 cos (r-|~3/2 _ ^-3/2) /y t 


but if ri ^r 2 , then there exists a sequence of t values 



with lim tn = 00 

rVk-m 11 


such that 

cos ((ri~3/2 _ r 2 ~3/2) /y t n ) = 0 


in fact, set 


f r 1 _ 3/2 _ r 2 ~ 3 / 2 \ /y t n = 7 + nfl , n = 0 , 1 , 2 ,... 


Thus 


I I Hi ( t n ) - D 2 tt n . 5 1 = /2, n = 1 , 2 ,... 

and therefore, if e < /2 is considered, the instability of the circular 
trajectories has been established. 
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3.0 LINEAR SYSTEMS 

Linear systems comprise the most important subset of the differential system 
equation (1), especially because the linearization of equation (1) about a 
particular solution is an important tool in the analysis of the originax 
system. Consider the linear system 


y = A(t)y (5) 

where A(t) is a continuous n by n matrix in the interval (.t 0 ,°°). The 
characteristic polynomial of equation (5) is the polynomial of degree 
n in A, 

| A ( t ) - \l| = 0 (6) 

where | | denotes the determinant of A(t) - XI. 

If A(t) is a constant matrix, then the constant coefficient system 

y = Ay (7) 

has a stable characteristic polynomial if every eigenvalue X has Re(A) < 0 
where Re denotes the real part of the possible complex number X. 

A very useful result for arbitrary linear systems is the next theorem, which is 
proved in reference 1. 

Theorem 1 : A necessary and sufficient condition for stability of every solution 

of the linear system equation (5) is that all solutions of equation (5) be 
bounded. 

A frequently occurring case for orbital mechanics is where the linear system 
equation (5) has multiple, purely imaginary eigenvalues X = +iB (usually 
occurring an even number of times), or if zero is a multiple eigenvalue. This 
case can be either stable or unstable depending_on the nature of the solutions. 
One procedure is to find a fundamental matrix £(t) of equation (5) and apply 
ireorem 1. A fundamental matrix is defined as an n by r. matrix such that. 

£(t) = A(t) £(t) 


and 


det £(t) i 0 , t 0 < t < 00 
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(i.e., it contains n linearly independent (vector) solutions of equation (5)). 
Example 2 (unstable). The system 


I = A Z» where A 


( 0 0-1 0 \ 
0 0 10 ] 
0 0 0 0 / 
0 0-1 0 / 


(E2.1) 


has the characteristic equation |a - Xl| = X 1 * = 0 and, hence, X = o is an 
eigenvalue of multiplicity 4. To determine the fundamental matrix for 
this system, eigenvectors £ = (ci ,02,03,04^ associated with X = 0 are 
sought. In this case, the system of four equations 


A£ = X£ 


has three linearly independent solutions 


£l = (1,0,0,0) T 
$2 = (0,1,0,0) T 


and 


£3 = (0, 0,0,1)? 

Consequently, £j.e^t = £^ are th^ee linearly independent solutions of equation 
(E2.1). To find a fourth, linearly independent solution of j = A^, £4, which 
satisfies (A - Xi) 2 £4 = 0, but does not satisfy (A - Xi )£4 = 0, is sought. 
Such a vector is £4 = (0,0, 1,0)?, and the corresponding solution is 

e^ t (£4 + t( A - XI) £4) = ( — t , t , 1 , — t ) T 

(This is easily verified by substitution into the aquation.) A fundamental 
matrix for this system is 
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t ) 


where 


e Xt 



0 

1 

0 

0 


0 -t 
0 t 
0 1 


1 -t, 


eXt = •) 


(E2.2) 


The fourth column of £(t) is a solution of equation (E2.1) and is unbounded 
because 


1 1 (-t,t, 1 ,-t) T | | = V3t 2 + 1 


Thus, Theorem 1 verifies the conclusion that equation (E2.1) is unstable. For 
further details on the generation of the solutions of linear systems with 
multiple eigenvalues, see reference 2, section 3.10. 

Example 3 (stable) . Consider the two oscillatory systems 


x, + u>i 2 xi = e-| Cx ! ,xt ,X2,x 2 )i 

X2 + a)2 2 X2 = E2 ( X 1 » X 1 » x 2 > x 2) ’ (E3.1) 


where oo q and k >2 are constants. Equation (E3-1) can be expressed as a sys- 
tem of first-order equations by the usual substitutions 


X 1 = PI 
Pi = -^1 2x 1 
x 2 = P2 
p 2 = -a >2 2x 2 


e 1 Cx i , pi ,X2,P2) 


€2( X 1 ,P1 , X 2»P2) 


I 

) 


CE3.2) 
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or in matrix form 



e = (0,e 1 ,0,e 2 ) T 

and A denotes the matrix in equation (E3-3). then system equation (E3.3) can 
be written. 


x = Ax + e 


(E3.H) 


For the present, £ is suppressed by setting it equal to zero. The resulting 
linear system 


Ax 


(E3.5) 


has characteristic equation 


A -XI = (A^ + lo-|^)(X^ + u) 2 ^) = 0 


and eigenvalues 


X = +iu) ■) 
X - +iu) 2 


If 4 . u> 2 , then it is well known (see Theorem 3) that system equation 

(E3.5) is stable; i.e., every solution of equation (E3.5) is stable. Consider 
the case where = u> 2 so that li w i are eigenvalues of multiplicity 2. 
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Proceeding as in equation (E2.1), linearly independent eigenvectors 
£ = (c-| ,c 2 ,C 2 ,C 4 ) t corresponding to X = iu) -j are sought, such that 

A£ = iu)^ (E3.6) 

Solving equation (E3-6) for £, it is found that £1 = ( 1 ,iWi ,0,0) T and 
£2 = (0 ,0 , 1 , ioo 1 ) T are linearly independent solutions of the algebraic system 
equation (E 3 . 6 ), and the corresponding solutions to equation (E3.5) are given by 

x 1 = e 1U)lt ( 1 ,iu) n ,0,0)T 

and 

X 2 = e iU)2t ( 0 , 0 , 1 , iu 1 ) T 
Similarly, with X = -iUi, solve 
A£ = -iu-|£ 
to find solutions 

X 3 = e iWlt ( 1 ,-iu-| , 0 , 0 ) T 
and 

X4 = e _iU3lt (0,0,1 ,-1^1 ) ^ 

A fundamental matrix is 

<Rt) = (x^t), X 2 (t), X3U), x 4 ( t ) ) 

Moreover, because each solution Xi^t), i = 1,2,3,** is bounded, it is concluded 
that the linear system equation (E3.5) is stable. 
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Example 4 (unstable) . Reference 5 considers the Kustaanheimo-Stiefel (KS) el- 
ement formulation of the two-body problem without perturbations. The linearized 
equations are 

» 

6y = A5y ( E4 . 1 ) 

where 



where E is the eccentric anomaly used as the independent variable and 


6^ = (60^, 6a 2 , 603, 604, 6B 1t S8 2 , <$83, <$84, 6u>, 6 t) t 


In this case, |A - X 1 1 = X 10 so that X = 0 is an eigenvalue of multiplicity 
10. Solutions to equation (E4.1) can be generated by first finding eigenvectors 
^ : (c;, c 2 , ..., 09 , cio)^ 3uch that 

A L = K 


which, in this case, reduces to 



0 


Thus, C9 = 0 (as - — /0) and nine linearly independent eigenvectors are 

800^ 

given by §.i = (0, ..., 0, 1,0, ... 0 ) T where 1 is in the i th position, and 
zero otherwise for i = 1, ..., 1C; i *9- The tenth eigenvector (denoted £9) 
can be found using the technique of the previous examples, and is 
£9 =(0, ..., 0, 1, 0) T . The corresponding solution to equation (E4.1) is 
given by 


6jr 9 = £ g 


E (A - XI) 



0 , 


3ME\ T 
80) V 
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The other nine linearly independent solutions to equation (E4.1) are given by 


^ = £i, i = 1, ..., 10; i i9 

Because 6yq is unbounded as a function of E, it is concluded that equation 
(E4.1) is unstable tin the sense of Lyapunov). 

The conclusions of the preceding examples could have been brought forth much 
faster by discussion of the rank of the matrix (A - XI), which indicates the 
number of linearly independent solutions of the linear system 


(A - XI) £ = 0 


however, the constructive approach seemed more appropriate. Because of the im- 
portance of the situation encountered in these examples, the following remark is 
stated. 

Remark . Suppose that in a real, linear system £ = ky_, for every eigenvalue X 
such that ReX = 0 and X has multiplicity m, the equation 


A L = xg. 


has m linearly independent solutions (g.i}i=i then i. = A^ will be stable 
or unstable in accordance with whether the remaining solutions (associated with 
the other eigenvalues) are bounded or unbounded. This is a consequence of the 
fact that the solutions ^ = e At Ej^, i = 1, . .., m are bounded. 

A second point to be aware of is that it suffices to study only the trivial, or 
zero, solution of linear systems, because by linearity 


Hl(t;t 0 ,Il) - Zo(t;t 0 ,£o) = n(t;t 0 ,£i - £ 0 ) 
is a solution, and conditions (S2) and (SI) reduce to 
I - ZoH I e. t > t Q 


whenever 


IZl -loll <5 
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Another implication of this second point is that constant coefficient linear sys- 
tems can be classified as stable, asymptotically stable, or unstable because 
whatever is true for the zero solution is true for all solutions. This is in 
contrast to nonlinear systems where specific solutions must be concentrated on; 
i.e., boundedness and stability are independent concepts, in general. 

Theorem 2 : If the characteristic polynomial (eq. (6)) of the linear system (eq. 

(7)) is stable, then every solution of the linear system (eq. (7)) is asymptot- 
ically stable. (See reference 1 for verification.) 

The word stable is unfortunately associated with the characteristic polynomial, 
because it implies more than that for the solutions. A final result in this di- 
rection is as follows: 

Theorem 3 : If the roots of the characteristic polynomial (eq. (6)) of the linear 

n 

system (eq. (7)) are 1 X i> ^ (not necessarily distinct), then every solution of 

the linear system is stable if (1) Re(X^) < 0 for all eigenvalues X^ with 
Re(X^) ^0, and (2) if there exist roots X^ with Re(X^) = 0, then X^ is a 
simple root (i.e., has multiplicity 1). 

Remark . Consider the constant coefficient system ^ = A^. Suppose A is unsta- 
ble because of the existence of an eigenvalue X with ReX >0. If there exists 
a nonsingular matrix M that diagonalizes A; i.e., 


M-UM = D = diag (X 1f ..., X n ) 


then it is well known that 


| A - Xl| = | D - Xl| 

and consequently the transformed system x = Dx is unstable where x = M -1 ^. 
This implies that if instability due to ReX >0 in a system i = Ay_ is 
eliminated, it will be necessary to use a change of independent variables or 
some transformation other than a similarity transform. 


3.1 PERTURBED LINEAR SYSTEMS 

Several of the preceding theorems were applicable to the constant coefficient 
linear system £ = A^, with the exception of the general nonautonumous linear 
system £ = A(t) Perturbations of asymptotically stable (or stable) con- 
stant coefficient linear systems can inherit stability provided the perturba- 
tions are small in some normal sense. Specifically, the linear systems under 
consideration are 
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1 = (A + C(t))jr 


( 8 ) 


where C(t) = (Cij(t))i, j = i is a continuous n by n matrix on (t 0 ,°°). In what 
follows , 


1 1 c ( t ) 1 1 = (SScfj(t)) 1/2 


or some other finite dimensional norm such as 



provided such maximum exists. 

Theorem 4 : If the characteristic polynomial of jr = is stable (i.e., ^ 

is asymptotically stable) and 


/ 1 1 C( t) 1 1 dt < » 

to 


(9) 


(i.e., the integral exists), then all solutions of equation (8) are asymptotically 
stable. (See reference 1 for verification.) 

Corollary 1 (to the proof of Theorem 4) . If there exists a sufficiently small 
constant fj > 0 such that 

| | C( t) | | < £, t included in the interval (t 0 ,°°) (10) 


then the conclusion of Theorem 4 is valid. The exact nature of £ is discussed 
in reference 1. Observe that inequality (eq. (10)) does not imply inequality 
( eq . ( 9 ) ) . 

Corollary 2 . If all the 3olutions of ^ = Ay are stable and condition 
(eq. (9)) holds, then all the solutions of the perturbed system (eq. (8)) 
are stable. (See reference 1 for verification.) 
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H.O PERTURBED NONLINEAR SYSTEMS 

Analysis of the stability properties of solutions of nonlinear systems can often 
be achieved by considering the nonlinear system 

i = 

as a nonlinear perturbation of a linear system; i.e., 

y. = A^ + P(t,^) (11) 


This situation occurs when attempting to analyze a particular solution of sys- 
tem equation (1) by linearizing that system about the particular solution. 
Because of the importance of chis process, details are provided below. 

Linearization . Let y 0 (t) be a given particular solution of the nonlinear sys- 
tem equation (1). Study the stability of y G (t). Let y = y(t) be any other 
solution of equation (1), then using a vector valued Taylor series expansion of 
f (expanded about y^t)), the result is 


f(t,y(t)) 


f(t,y 0 (t)) + 


df 

— ( t ,y 0 ( t) ) 

- d Z 


ty “ Zo(t^ + 



((y - y 0 )-Z)' 


m=2 


m ! 


f(t ,y 0 ( t) ) 


( 12 ) 


/ 3 3 \ T 

where V = - — , - — , fU.y,-,) = ( f i ( t ,y 0 ) , .... fn^.Zo^ 1 and 

\ 3 yi 3 y n / 

((y - y D ) -Vjm is a symbolic operator in which 


(i - I„> l 


03) 


is performed and then raised algebraically to the m 111 power. At this point, it 
is applied to the components of f and evaluated at (t,y 0 (t)). Next, observe 
that 


f(t,y(t)) = y(t), f( t ,y Q ( t) ) = y 0 (t) 
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which implies that equation (12) can be written as 


51 (I - *>' = 


df 

- (t,J[o) 


‘•Z - 


4E 


<<z - Zo>*Z> m 


^m=2 


m! 


I(t,£ 0 ) 


Set 


u = x ( t ) - £ 0 (t) 


(14) 


in the last equation to obtain 


u 


3f 

— (t,£o(t)) 

3 Z 



f(t,£ 0 (t)) 


(15) 


or, 


u = Au + £(t,u) 


(16) 


where 


A 


3 f 

— (t,£ 0 (t)) 
3 Z 


the Jacobian of f at £ 0 (t) 


(17) 


and 


£(t,u) 




l(t,jr 0 (t)) 


( 18 ) 


The linear part of equation (16), which is obtained by suppressing £(t,M), 
is called the linearized system associated with equation (1), and is often 
written 
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$Z = 



«z 

U,y 0 (t)) 


(19) 


In some instances, a linearization about initial conditions (t 0 ,£ 0 ) is useful, 
and in this case the result is 


6* = 


"3 f " 

. 3 Z. 


iz 

( t 0 ,^ 0 ( t) ) 


( 20 ) 


Some of the most important cases of linearization occur for autonomous systems 


z = Hz) 


( 21 ) 


The preceding formulas are modified accordingly. The following examples will 
illustrate most of the ideas of this section. 

Example . Consider the nonlinear equation from reference 3- 


y = i - y 2 , y(t 0 ) = y Q (E5.D 

The real solutions of equation (E5.1) are given by 
y = tanh (t - t 0 ♦ k) 

where 


k = arc tanh (y Q ), -1 < y 0 < 1 (E5.2) 


It is easy to show directly that the solution y 0 = -1 is unstable as t ■* OD , 
whereas yi = 1 is stable. The objective is to relate equations (i), (16), and 
(19) for this example. Here, 


f(t,y) = 1 - y 2 
3f 

— (t,y) = -2y 

dy 
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and if an expansion occurs about y Q = -1, 


A = 


3f 

3 y. 


= 2 


(E5.3) 


yo 


Moreover, regarding p(t,u), the following occurs 


(u-V)® j-u 2 , if m = 2 

( 0, if m > 3 


m! 


(E5.4) 


Consequently, equal: on (16) equals 


u = 2u - u* 


(E5.5) 


where 

u = y-yo(t)=y + 1 (E5.6) 

The linearized equation (19) is 

6y = 26 y (E5.7) 

Assuming 0 < u < 2 (for convenience), then equation (E5.5) can be solved, and 
2ce 2 ^ 

u(t) = , c > 0 ( E5 . 8 ) 

1 + ce 2t 


As noted by l'Hopital's rule that u(t) *2 as t 00 , thus because u = y ♦ 1 , 
this implies that y(t) 1 as t ♦ 00 . Thus, every solution starting near 
y 0 ( t) = -1, but greater than -1, approac 1 ;s asymptotically y| : 1. Hence, 
y Q = -1 is unstable. Clearly, the linearized system equation (E5.7) has 
solutions 


6y 


= ke 2t 


(E5.9) 
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Figures 3. 4 , and 5 illustrate various solutions of systems equations (E5.1), 
(E5.5), and (E5.7), respectively. Figure 4 is the same as figure 3 but 
translated one unit upwards (in the y direction). Figure 5 (illustrating 
the linearized system) only fits the nonlinear system initially. In this 
case, the nonlinear part - u ? results in the solutions being asymptotic to 2. 

A similar expansion could be made about yi = 1. For this example, the instabil- 
ity of the linearized system correctly predicted the instability of the solution 
y 0 = -1 of the nonlinear system. Two important remarks are in order. 

Remark 1 . The linearized system equation (19) will not. always correctly predict 
the stability behavior of the nonlinear system equation (1). It is possible 
for equation (19) to be stable but for the nonlinear system equation (1) to 
have a solution y 0 (t), which is stable or unstable depending on the nature 
of £(t,u) in equation (16). Essentially, a stable situation in which there 
are multiple eigenvalues of linearized systems of the form X = +i6 can become 
stable or unstable depending on the original nonlinear system. This situation 
is illustrated with a classic example in reference 2. (See equation (E6.1)). 





y = 1 - y 2 


Figure 3.- Solutions 

. . . 2 
of y ■ 1 - y 


u = 2u - u^ 
where 

u=y-y 0 =y*i 

Figure 4.- Solutions 
of u - 2u - u" 


<5y = 2<$y 

superimposed over 
solutions of equation 
(E5.5) 

Figure 5.- Solutions 
of (y “ 26y 
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Remark 2 . When studying autonomous nonlinear systems and the stability of cer- 
tain solutions, many times a change of dependent variables is made that essen- 
tially translates the system to the origin. This same result is achieved by 
linearization. 

For example, the nonlinear system 


y = 1 - y 2 

becomes 


u = 2u 


2 


under the change of variables u = y + 1 . 
Example 6 . Consider the nonlinear system 


yA / yz - y i (y i 2 + y2 2 > N 
,y 2 / \-y i - y2<yi 2 + V2 2 h 


(E6.1) 


which has ^(t) = ( y *; ( t ) , y 2 ( t ) ) ^* =0 as a solution. The linerrized system is 


^yr 

V 6 y 2 ) 


0 r 

k-1 0> 


^yi' 

< & Y2/ 


which has eigenvalues A r +i. Thus, the linearized system is stable 
observe that 


(26.2) 
Next , 


— (yi 2 + y 2 2 ) = -2 ( y i 2 + y 2 2 ) 2 
dt 


(E6.3) 


where equation (E6.1) is used to evaluate on the left 2 ( y -) y -j + y2y2) • Setting 


r(t) = yi 2 + y 2 2 


equation (E 6 . 3 ) equals 
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which integrals to give 


1 

= 2t + c (E6.4) 

r(t) 


This last equation implies 


yi 2 (t) + Y 2 2 (t) = 

1 + 2ct 


(E6.5) 


where 


n = y-| 2 (o) + y 2 2 (o) > 0 


(pro. 'ded the initial value problem is nontrivial). Recalling that 
Hy.il = Vy-[ 2 (t) + y 2 2 (t) , equation (E&.5) implies that 


lim ||j_|| = 0 


where £ = y.(t;t 0 ,£ 0 ) and y_ 0 = (vi (o) ,y 2 (o) ) T . Thus, the result is that 
0 is an asymptotically stable solution of equation (E6.1). Because, in 
addition, 


I lz(t) 1 1 < 1 |y.(o) 1 1 , t > 0 . 


Example 7 On 



the other hand, 

( >’2 + yi^yi 2 
-yi - y2(yi 2 


isider the system 

y2 2 )\ 

y2 2 )/ 


(E7.1) 
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The linearized system associated with equation (E7.1) is exactly the same as 
equation (E6.2), which, in turn, is stable. In this case 


d 

dt 


(yi 2 + Y2 2 ) = 2 (yi 2 ♦ Y 2 2 ) 2 


(E7.2) 


which implies 


2 2 

y i + y2 = 



(E7.3) 


where 


c = yi 2 (o) + y 2 2 (o) 

1 

The right side of equation (E7.3) has a vertical asymptote at t = — > 0; 

2c 

i.e., the solution has finite escape time. Thus, 


lim I ly(t;t 0 ,y_ o ) 1 1 = 00 
trM/2c 


regardless of c. This implies that equation (E7.1) is unstable at the origin. 

The examples presented above indicate that stability of the linearized system 
will not necessarily be inherited by the nonlinear system. The following theo- 
rems indicate more positive results for 

Z = Ajr + £(t,z> 

The following hypothesis regarding the system equation (11) is required: 

a. £(t,z) is continuous in T, and 

li£(t, Z )|| 

b. lim = 0, uniform in t. 

I III l-o llzll 

Condition (b) implies that £(t,o) = o so z = H is a solution of equation 

( 11 ). 
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Theorem 5 . Suppose |a - Xl| is stable (i.e., the linear system is 
asymptotically stable) and £(t,^) satisfies (a) and (b), then the solution 
^ = o of system equation (11) is also asymptotically stable. (See reference 
1 for verification.) 

A similar version for autonomous systems 


£ = Ay_ + £(£) 


( 22 ) 


appears in reference 2, where the norm 



= max | yjj 
i 


is used. 

Theorem 6 . Suppose the vector valued function 


£(£) 

TiliT 


(23) 


is a continuous function of (i.e., the components y^ of are continuous) 
and vanishes for £ = 0. Then 

a. The solution = jD of equation (22) is asymptotically stable if the 
linearized system <$£ = A 6y is asymptotically stable. 

b. The solution £ = 0 of equation (22) is unstable if there exists at least 
one eigenvalue X of the characteristic polynomial |A - Xl | such 

that ReX > 0. 

c. There is no conclusion regarding the stability of = 0 if every 
eigenvalue X of the characteristic polynomial |a - X I | has ReX < 0, 
and there is at least one eigenvalue with ReX = 0; i.e., there exists a 
purely imaginary eigenvalue of A. 

Example 8 . For the final example, consider the one-dimensional perturbed 
harmonic oscillator, 


y 


+ £ ( 1 - y 2 )y + y = 0 


( E8 . 1 ) 
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where e is a small, positive parameter. This equation is equivalent to the 
system 


( 




(E8.2) 


where y-| = y» Y2 = y* The linearized system is 



(E8.3) 


which has eigenvalues (distinct for 0 < e < 2) 


X = 



2 


Because ReX < 0, the linearized system is asymptotically stable to (0,0)^. 
Moreover , 


£( v) 

lim = 0 (E8.4) 

z*o llzll 

where £(y.) = C 0 , ery -j ^y 2 ) T - Thus, by either of the preceding theorems, (0,0)^ 
is an asymptotically stable solution of equation (E8.2). 

With regard to equation (E8.4), it should be noted that it suffices to consider 


e yi^y2 


llzIN V yi 2 «■ y2 2 
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but 


and 


eyi 2 y2 


eyi 2 y2 y/y i 2 + y2 2 £ y2 \/yi 2 + y2 2 


Vyi 2 ♦ y ? 2 


2 2 

yi ♦ yi 


w( y -l 

\y2> 


eyg Vyi 2 y2 2 
i + (yi/y2 ) 2 


= < e |y 2 | Vyi 2 + y 2 2 


(E8.5) 


The inequality (eq. (E8.5)) implies the conclusion (eq. (E8.4)) 


5.0 CONCLUSIONS 

The study of the qualitative behavior of solutions of linear and nonlinear sys- 
tems of differential equations is of major importance and has implications for 
the accuracy of numerical solutions. It is important for scientists and engi- 
neers involved in the analysis (numerical and analytical) and the solution of 
large scale systems to clearly understand the concept of stability both from the 
numerical and mathematical viewpoint. For various reasons, most publications 
emphasizing the mathematical and numerical stability theory of linear systems 
tend to focus on small systems; i.e., one or two first-order differential equa- 
tions. This does not necessarily enhance understanding of the larger systems. 

In orbital mechanics, systems of 6 to 10 equations are more prevalent, and they 
frequently have high multiplicity eigenvalues, which are either zero or purely 
imaginary. In this report, the emphasis is on mathematical (Lyapunov) stability 
analyses for these higher dimension systems with eigenvalues of the form ifj 
where i 2 = -1 and 0 are real, nonnegative numbers. Eight examples are 
presented from orbital mechanics, mechanical, and electrical systems, etc., 
which illustrate the techniques and theory of Lyapunov stability. 

A basic tool in the study of nonlinear systems of equations is to linearize the 
system and analyze the resulting linear system. Unfortunately, she stability 
properties of the linear system are not necessarily inherited by the original 
nonlinear system, and this is illustrated by example. Of more concern is that 
if the linearized system is only (marginally) stable (i.e., it is stable but has 
repeated eigenvalues of the form i0), or (marginally) unstable, then the 
nonlinear system may go either way. This is precisely the situation often 
encountered in orbital mechanics. This report attempts to emphasize understand- 
ing of Lyapunov stability by illustrating techniques and giving examples that 
occur in practice, and to stress the type of system for which analysis is more 
difficult. 
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